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Calculations  for  Uniaxial  Stress  and  Strain  Cycling 
in  Plasticity 


by 


P.  M.  Naghdi  and  D.  J.  Nikkei,  Jr. 
Department  of  Mechanical  Engineering 
University  of  California 
Berkeley,  CA  94720 


Abstract .  Within  the  framework  of  an  existing  purely  mechanical,  rate-type 
theory  of  plasticity  detailed  calculations  are  presented  for  certain  types 
of  material  response  during  stress  and  strain  cycling  in  a  uniaxial 
homogeneous  deformation.  These  features  pertain  specifically  to  material 
response  in  stress  cycling  between  fixed  values  of  stress  in  tension  and 
compression  (not  necessarily  equal  in  magnitude)  resulting  in  ratcheting 
of  strain,  and  a  type  of  saturation  hardening  caused  by  strain  cycling 
between  any  two  fixed  values  of  strain  when  the  mean  value  of  stress  (in 
tension  and  compression)  tends  to  zero. 


1. 


Introduction 


Starting  from  a  general  theory  of  plasticity  developed  in  stress  space 
by  Green  and  Naghdi  [1,2]  and  its  formulation  relative  to  strain  space  by 
Naghdi  and  Trapp  [3],  Caulk  and  Naghdi  [4]  have  discussed  hardening  response 
and  in  particular  have  compared  the  predictions  of  the  theory  with  two  sets 
of  experimental  results.  In  [4],  with  the  limitation  to  small  deformation, 
systematic  specialization  of  the  constitutive  equations  of  the  general  theory 
were  specifically  applied  to  saturation  hardening  and  other  types  of  stress 
or  strain  cycling  were  not  considered.  Additional  developments  of  saturation 
hardening  in  the  context  of  strain-space  formulation  were  discussed  by  Casey 
and  Naghdi  [5];  and,  for  the  class  of  materials  for  which  detailed  compari¬ 
sons  with  experiments  were  made  in  [4],  the  results  for  saturation  hardening 
in  [4]  and  [5]  are  equivalent. 

More  recently,  the  stress  and  strain  cycling  has  also  been  discussed  by 
Drucker  and  Palgen  [6],  who  employed  an  approach  to  plasticity  different  from 
that  of  strain-space  formulation  utilized  in  [4,5].  In  the  course  of  their 
development,  Drucker  and  Palgen  [6]  have  also  discussed  the  material  response 
appropriate  for  unsynmetric  plastic  cycles  of  stress  or  of  strain.  While  it 
may  not  be  obvious  at  first  sight,  one  would  expect  that  a  suitable  specializa¬ 
tion  of  the  general  theory  should  be  capable  of  accommodating  the  type  of 
unsyranetric  cycles  mentioned  and  discussed  by  Drucker  and  Palgen  [6].  The 
purpose  of  this  paper  is  to  show  that  this  is  the  case.  In  particular,  in 
what  follows  we  demonstrate  the  predictive  capability  of  the  theory  in  [1,2,3] 
for  two  types  of  uniaxial  material  response,  namely 

1)  ratcheting  of  strain  caused  by  stress  cycling  between  fixed 
values  of  stress  in  tension  and  compression  (not  necessarily 
equal  in  magnitude)  similar  to  that  shown  in  Fig.  1  and 


2)  progressive  relaxation  to  zero  of  the  mean  stress  caused  by 
strain  cycling  between  any  two  fixed  values  of  strain  similar 
to  that  shown  in  Fig.  2. 

The  above  two  types  of  responses  in  the  light  of  experimental  results 

were  described  and  discussed  by  Drucker  and  Palgen  [6,  p.  480].  By  way  of 

background  and  for  definiteness,  some  features  of  these  two  types  of  responses 

are  recalled  here  with  reference  to  a  uniaxial  homogeneous  deformation  in  a 

tension-compression  test.  In  the  context  of  small  deformation,  let  e  and  s 

stand,  respectively,  for  the  component  e^  of  the  strain  tensor  and  the 

* 

component  s^  of  the  stress  tensor  .  Then,  a  stress  cycle  such  as 

AjBjCjDjE in  Fig.  1  consists  of  loading  in  one  direction  from  some  chosen 

value  of  stress,  say  at  ,  to  a  limiting  value  of  stress  at  Cj  followed  by 

reversing  the  direction  of  loading  until  another  limiting  value  of  stress 

(not  necessarily  equal  to  that  at  is  reached  at  Fj  and  the  direction 

of  loading  is  again  reversed  until  the  stress  reaches  a  value  at  A^  equal 

to  that  of  the  initial  choice  at  Aj.  More  generally,  again  with  reference 

to  Fig.  1,  during  the  n1*1  cycle  (n  =  1,2,3, .. .)  the  response  which  begins  at 

any  point  such  as  A  ,B  ,...,  or  F  and  ends  at  a  corresponding  value  of 
n  n  n 

stress  at  A  ,,B  or  F  ,  constitutes  a  stress  cycle.  Consider 

n+1  n+1  n+1 

now  a  specimen  cyclically  loaded  into  the  plastic  range  between  stress 

§ 

limits  whose  average  is  nonzero  i.e,,  stress  limits  which  are  unsymmetric 
about  the  strain  axis  in  the  s-e  plane.  During  a  typical  stress  cycle  of 
this  kind  (e.g.,  A^B^CjDjH^F^A,  in  Fig.  1)  the  strain  at  the  beginning  and 
end  of  the  cycle  (and  hence  also  the  minimum  and  maximum  values  of  strain 
attained  during  the  cycle)  vary  from  one  cycle  to  the  next.  Evidently, 

* 

A  distinction  between  the  Cauchy  stress  tensor  or  either  of  the  Piola- 
Kirchhoff  stresses  is  unnecessary  when  only  a  small  deformation  is  being 
considered. 

§ 

This  feature  is  referred  to  in  [6]  as  "unsymmetric  cycles  of  stress  in 
the  plastic  range." 


2. 


experimental  observations  (see,  for  example  Dolan  [7,  Fig.  14])  indicate 

* 

that  if  the  average  of  the  stress  limits  is  positive,  both  the  minimum 
and  maximum  strains  attained  in  a  given  cycle  can  progressively  increase 
from  one  cycle  to  the  next.  Similarly,  it  may  be  inferred  from  the  dis¬ 
cussion  in  Sandor's  book  [8,  pp.  97-100]  that  if  the  average  of  the  stress 
limits  is  negative,  then  both  the  minimum  and  maximum  values  of  strain 
attained  in  a  given  cycle  can  progressively  decrease  from  one  cycle  to  the 
next.  Thus,  while  the  stress  limits  remain  fixed  for  repeated  cyclic 

loading,  the  strain  limits  of  the  cycles  can  be  progressively  displaced 

** 

along  the  strain  axis  from  one  cycle  to  the  next  .  A  well-known  phenomenon 

of  this  type,  called  ratcheting  of  strain,  refers  to  progressive  increase 

in  plastic  strain  during  cycling  between  any  two  fixed  values  of  stress 

levels.  Evidently  such  a  progressive  increase  in  plastic  strain  may  take 

place  irrespective  of  the  starting  point  in  the  cycle.  Indeed,  as  is  clear 

from  Fig.  1,  the  strains  at  A  ,,B  F  ,  all  have  been  increased  along 

6  n+1’  n+1  n+1  & 

the  e-axis  from  the  values  of  strain  at  A  ,B  ,...,F  . 

n  n  n 

For  cycling  in  the  plastic  range  between  fixed  strain  limits,  on  the 
other  hand,  for  most  metals  the  extent  of  the  elastic  range  during  each 
successive  strain  cycle  increases  (or  decreases)  until  a  limiting  periodic 
response  is  reached  in  which  the  stress-strain  curve  of  each  succeeding 
cycle  is  the  same  (see  Fig.  3).  This  phenomenon,  which  was  previously  dis¬ 
cussed  in  [4],  is  sometimes  called  saturation  hardening  and  a  material  is 
said  to  saturate  when  it  reaches  this  limiting  behavior.  Moreover,  if  the 
★ 

The  terminology  of  the  average  of  the  stress  level  is  referred  to  in  [6] 
as  the  "mean  stress"  and  "mean  stress  in  the  cycle." 

The  strain  limits  of  the  cycles  will  be  displaced  in  the  positive  strain 
direction  if  the  average  of  the  stress  limits  is  positive  (as  depicted  in 
Fig.  1)  and  in  the  negative  strain  direction  if  the  average  of  the  strain 
limits  is  negative  (see  [8,  Fig.  98]). 


1 


specimen  is  cycled  between  fixed  strain  limits  whose  average  is  nonzero,  it 

§ 

has  been  observed  by  Morrow  and  Sinclair  [9]  and  by  Dolan  [7,  Fig.  15] 
that  for  a  certain  class  of  materials  the  average  of  the  limiting  stresses 
in  the  cycle  will  progressively  relax  and  tend  to  zero.  For  example,  with 
reference  to  Fig.  2,  the  average  value  of  the  stresses  at  Cn  and  Fn  (n  >  3) 
tends  to  zero  for  repeated  cyclic  loading  between  fixed  strain  limits. 

For  completeness,  we  also  briefly  indicate  in  Fig.  3  the  prediction  of  the 
theory  (by  means  of  the  procedure  of  the  present  paper)  for  the  case  of 
saturation  hardening,  and  examine  its  comparison  with  that  in  [4]. 


I 


See  also  [6, 


Fig.  2]. 


4. 


2. 


Basic  Equations 

Using  a  strain-space  formulation  of  plasticity,  we  summarize  here  the 
basic  equations  of  a  purely  mechanical  theory  contained  in  the  papers  of 
Green  and  Naghdi  [1,21,  Naghdi  and  Trapp  [3]  and  Casey  and  Naghdi  [5).  In 
addition  to  the  strain  tensor  e^>  we  assume  the  existence  of  a  symmetric 
second  order  plastic  strain  ej^  and  a  scalar  k  representing  a  measure  of 
work -hardening.  It  is  assumed  that  the  stress  tensor  svt  is  given  by  a 
constitutive  equation  of  the  form 


SKL  =  SKL(-U-1 


U  *  (eMN’eMN,IC) 


(1) 


and  that,  for  fixed  values  of  eJL.  and  k,  (1)  can  be  inverted  in  the  form 

lllNi  I 


eKL  ~  eKL^  *  V  *  fsMN,eMN,IC^  ' 


(2) 


The  response  functions  s^  and  e^  in  (1)  and  (2)  are  taken  to  be  smooth. 
We  admit  the  existence  of  a  smooth  yield  (or  loading)  function  g(U)  in 
strain  space  such  that  for  fixed  values  of  the  equation  g(U)  =0 

represents  a  closed  orientable  hypersurface  of  dimension  five  enclosing  a 
region  of  strain  space.  With  the  use  of  (2)^,  we  obtain  a  corresponding 
function  f(l/)  through  the  formula 


g(U)  =  g(eMN(!T),ePN,K:)  =  f(l/)  .  (3) 

Because  of  the  assumed  smoothness  of  (2)  ^ ,  for  fixed  values  of  (e^»K) 
the  equation  f (10  =0  represents  a  hypersurface  in  stress  space  having  the 
same  geometrical  properties  as  the  hypersurface  in  strain  space. 

The  loading  criteria  of  strain- space  formulation  will  be  regarded  as 
primary.  Also,  for  our  present  purpose,  it  is  sufficient  to  consider  only 
a  special  case  of  the  stress  response  in  which  s^  in  (1)^  is  independent 


5. 


of  <  and  its  dependence  on  the  variables  occurs  only  through  the 

difference  e._, -ef...  Then,  the  constitutive  equations  for  the  rate  of 
MN  MN 

plastic  strain  ePL  and  rate  of  work-hardening  parameter  k  during  loading 

A  * 

(i.e.,  when  g  =0  and  g  >0)  may  be  expressed  as  (see  the  development  between 
Eqs .  (36) - (42)  in  [5]): 


= 


3f 


'KL  r+A  3s 


KL 


<:  =  C  ep 
K  LKL  KL 


(4) 


where 


8 "  ^ eKL  ' 

.  r  •  -  (-2X-.  M-c,  ,)  ,  r  ■*  a  >o 


3SKL  aeKL 


3k  KL 


(5) 

(6) 


and  where  the  response  coefficients  C  are  functions  of  the  variables  (1)-,. 

It  should  be  noted  that  in  obtaining  the  expression  (4)^  a  condition  resulting 
from  a  work  inequality  proposed  in  [10]  has  been  also  imposed  on  the  con¬ 
stitutive  equation  for  plastic  strain  rate. 

The  characterization  of  strain-hardening  behavior  by  means  of  the  rate- 
independent  quotient  f/g,  with  f  =  (Sf/Ss^s^,  is  discussed  in  [5].  More 
recently,  a  scalar  quantity  $  is  defined  in  [11]  in  such  a  way  that  does  not 
presuppose  a  condition  of  loading  and  is  such  that  during  loading  it  has 

A  A 

the  same  value  as  f/g.  For  our  present  purpose,  the  scalar  function  4>  can 
be  defined  as  [12,  Eq.  (4.13)] 


$  = 


r+A 


(7) 


and  a  particle  of  an  elastic-plastic  material  is  then  said  to  exhibit  [11,12] 


The  scalar  g  is  defined  by  Eq.  (5)  below. 


6. 


hardening  behavior  if  and  only  if  $  >0  ,  (a) 

softening  behavior  if  and  only  if  4>  <0  ,  (b)  (8) 

perfectly  plastic  behavior  if  and  only  if  $  =  0  ,  (c) 

Moreover,  the  strain-hardening  response  of  a  material  is  said  to  saturate 
to  one  of  the  values  {K^,Ks,0}  and  exhibits  either  (a)  saturation 
hardening,  (b)  saturation  softening  and  (c)  saturation  to  perfectly 
plastic  behavior,  in  that  order,  if 


lim  Q 
t-»° 


(a) 

(b) 
Cc) 


(9) 


7. 


3. 


Special  Constitutive  Response  Functions 

It  is  convenient  to  prescribe  the  various  constitutive  response  functions 
in  terms  of  the  components  of  deviatoric  stress  and  deviatoric  strain  ykl» 
namely 

(10) 

with  similar  definition  for  tne  components  of  deviatoric  plastii  rain 

4,  where  s  and  e  are  the  mean  normal  stress  and  mean  normal  si  rs. 

respectively.  In  the  context  of  small  deformation  and  for  mateT  ,  which 
are  isotropic  in  the  reference  state,  we  specify  the  stress  response  by 
generalized  Hooke's  law  and  this  can  be  written  in  the  form 

TKL  =  2p(YKl  _yKL5  ’  s=3k(e-eP)  ,  (11) 

where  U  is  the  shear  modulus  of  elasticity  and  k  is  the  bulk  modulus. 

In  regard  to  the  characterization  of  material  response  to  stress  and 
strain  cycling  of  the  types  referred  to  under  1)  and  2)  in  Section  1,  we 
first  observe  that  it  is  the  translation  of  the  yield  surface  (whose 
interior  represents  the  elastic  domain),  and  not  its  size,  which  is 
particularly  significant  here.  For  the  response  described  under  1)  in 
Section  1,  while  the  plastic  strain  progressively  increases  from  one  cycle 
to  the  next,  the  overall  translation  of  the  yield  surface  in  stress  space 
is  not  significantly  different  between  the  later  and  earlier  cycles.  For 
the  response  described  under  2)  in  Section  1,  on  the  other  hand,  the  amount 
of  translation  of  the  yield  surface  in  stress  space  becomes  progressively 
smaller  as  the  material  is  cycled  between  fixed  strain  limits.  In  the  yield 
function  utilized  previously  by  Caulk  and  Naghdi  [4,  Eq.  (40)],  the  quantity 
representing  translation  of  the  yield  surface  in  stress  space  can  be  regarded 


TKL  =  SKL  "  S  6KL  ’  TKK  =  ° 


YKL  =  eKL  "®  6KL  ’  YKK  =  ° 


as  a  constant  represented  by  h  a  times  the  plastic  strain.  Hence,  at  a  given 

value  of  plastic  strain  the  translation  of  the  yield  surface  is  determined 

independently  of  the  loading  history.  By  contrast,  in  our  present  discussion 

it  is  important  to  account  for  the  effect  of  loading  history  pertinent  to  the 

translation  of  the  yield  surface.  A  simple  way  to  achieve  this  is  to 

generalize  the  restricted  equations  of  Caulk  and  Naghdi  [4,  Eqs.  (40),  (56) 

and  (70)  ]  by  replacing  the  constant  coefficient  a  with  a  scalar  function  at 

of  the  work-hardening  parameter  <  which  is  0(1).  In  particular  we  take  a 

to  be  linear  in  k,  with  an  initial  value  a  when  <  =<  and  a  value  a  at 

oos 

saturation  when  k  =<s-  Thus,  we  specify  the  loading  function  in  the  form 


f  =  (T  -  “00  yP  )  ( X  -  yP  )  -  < 

1  KL  2  YKLnTKL  2  YKl/ 


[Y, 


(1  +  ^ykl][yki  ‘(1  +  ^IT)ykl]  =  g  ’ 


It  should  be  noted  that  for  the  special  case  in  which  a  ~a  =a,  a(<)  becomes 
simply  a  and  the  loading  functions  (12)  reduce  to  those  used  by  Caulk  and 
Naghdi  [4,  Eqs.  (40)^,  (56),  (70)^.  In  (17'',  h  may  be  interpreted 

as  the  center  of  the  yield  surface  in  the  deviatoric  stress  space. 

Similarly,  the  quantity  (1  +  a(<))y^  may  be  interpreted  as  the  center 

of  the  yield  surface  in  deviatoric  strain  space.  The  size  of  the  loading 
surfaces  (which  represents  the  extent  of  the  elastic  region)  depends  on  the 
work-hardening  parameter  k. 

Supplementary  to  (12)  and  (13),  we  also  need  to  specify  the  work¬ 
hardening  coefficients  in  (4).  As  noted  earlier,  a  significant  feature 


9. 


of  the  material  response  described  under  1}  and  2)  in  Section  1  is  the  trans¬ 
lation  of  the  yield  surface,  and  not  its  size.  However,  with  reference  to 
the  phenomenon  of  saturation  hardening  examined  in  [4],  it  is  the  size  of 
the  yield  surface  which  is  of  particular  interest.  In  order  to  account 

also  for  saturation  hardening,  it  is  sufficient  to  specify  in  the  same 

* 

form  as  that  in  [4,  Eq.  (56) ^],  by 


K-K 

C  =  ( _ - 

kL  V  -K 


*  (Vkl 


+  62YKL> 


-K  )  ^TKL  +nYKL)  * 
s  o 


(14) 


where  and  are  constants  and  where  with  the  use  of  (11)  B  and  n 

are  given  by  B^  =  2p8  and  8.,  =n  -2p6. 

With  the  constitutive  equations  (11)-(14),  it  follows  from  (6)  that  on 
the  yield  surface  the  scalar  functions  A  and  T  assume  the  values 


A  =  8vik  >0  ( f  =  g  =  0)  , 

A  A 

2coc  +  2  [1  +  (K  _K  )  (tkl  -  J  YKL)YKL^®TMN  +  nYMhP(TMN  *  2  YMN5 


a  -a 


(f  =  g  =  0)  , 


(15) 


(16) 


where  the  quantities  B  and  n  are  given  by 

/v  <  -<  ^  K-K 

8  =  B(k)  =  (^-~-)B  ,  n  =  nOO  =  (j-— )n  (17) 

os  OS 

subject  to  the  restrictions 


B(ks)  =  n(<s)  =  0  .  (18) 

In  view  of  (7),  (15)  and  (6)j,  it  is  easily  seen  that  the  strain-hardening 
★ 

The  hardening  response  coefficient  in  [4]  is  denoted  by  instead  of 
C„.  of  the  present  paper.  ‘  " 


10. 


is  characterized  by 


$  =  i 


8 

'  T 


11. 


w 


4.  Uniaxial  Cyclic  Loading 

Consider  now  a  homogeneous  extensional  deformation  sustained  by  a 
uniaxial  stress  s^ssft)  say.  Adopting  the  notations  e^  =  e,  ePj  =  ep, 
as  in  [4,5]  in  matrix  notation  we  write 


lit,,  II  4  »#>„  II  •  Ih-E,  II  ■  \  «p  llbKLll  , 


KL"  _3  ""  KL"  ’  11  KL"  2 

where  the  constant  matrix  ||bKLll  defined  by 


M  ■ 


2  0  0 
0-10 
0  0-1 


(20) 


(21) 


is  introduced  for  convenience.  For  the  uniaxial  homogeneous  deformation  under 
discussion,  each  of  the  two  variables  e  and  s  is  represented  as  a  line  in 
the  s-e  plane  (corresponding  to  the  six 'dimensional  strain  space  and  stress 
space) .  The  yield  surfaces  are  now  points  on  these  lines  and  the  elastic 
range  (corresponding  to  the  interior  regions  of  the  yield  surfaces)  is  a 
line  segment  between  these  points.  The  midpoint  of  the  line  segment  cor¬ 
responds  to  the  center  of  the  yield  surface  and  the  length  of  this  line 
segment  corresponds  to  the  size  of  the  yield  surface  whose  interior 
represents  the  extent  of  the  elastic  region. 

The  loading  function  (12)  now  assumes  the  simplified  form 

f  =  |  (s  -  |  a(«)eP)2  -  k  (22) 

* 

and  the  constitutive  equations  (11)  and  (4)j  2  with  (14)  reduce  to 

s  =  E(e  - eP)  (23) 

* 

We  have  used  the  form  (14)2  rather  that  (14)j  simply  to  make  the  form 
of  the  resulting  equations  more  compact. 
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and 


ep  = 


4  3  p  • 

3  E(s  -j  ae^)e 


(4  E  ♦a)  (s  aep)  ♦  [1  ♦  (-°-_~0  (s  otep)ep](6s  nep) 

os 


(24) 


K  =  (B(<)s  ♦  J  n(*)ep)ep  , 


(25) 


where  E  is  Youig's  modulus.  In  the  uniaxial  case  under  discussion,  A  is 

still  given  by  (15)  but  on  the  yield  surface  T  and  $  assume  the  values 

a  -a 
_o _ s 

c  -< 
o  s 

and 

r. 


F i  ■  2a<  +  i  [l  ♦  (-—."") (s  -  j  aep)ep] (6s  ♦  |  nep)(s  -  |  aep) 


$  = 


1 

ryTT  • 


(26) 


(27) 


Then,  from  (23) -(25),  we  have 

a  A  «  A 

dK  j(Bs-fnep) 

=  - p- -  . 

/S  A 

For  the  special  case  in  which  the  coefficients  a0  >“5*0,6  satisfy  the 
conditions 


(28) 


aQ  -  as  =  a  (say)  ,  2n  ♦a6  =  0  , 


(29) 


after  adopting  (29) 1  ^  and  making  use  of  the  condition  f  *0  with  f  given  by 
(22) ,  Tj  given  by  (26)  reduces  to 

T10  =  2(a+B)<  (30) 


and  (28)  becomes 


ck 

ds 


±26  (<) 

- T - 

a  ♦  8(<) 


(31) 
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where  in  writing  the  right-hand  side  of  (31)  again  use  has  been  made  of  the 
condition  f  *0  with  f  given  by  (22)  and  the  plus  and  minus  signs  correspond 
to  tensile  and  compressive  values  of  the  stress. 

In  view  of  (6)3  and  (7),  in  the  present  discussion  T  can  also  be  used 
to  characterize  the  strain -hardening  behavior.  Keeping  this  in  mind  and  the 
fact  that  k  >0,  it  follows  from  (30)  and  definitions  in  (8)  that 


a  +S(k) 


>  0  for  hardening  , 

<  0  for  softening  , 

=  0  for  perfectly  plastic 


(a) 

(b) 

(c) 


(32) 


The  results  (32)  hold  only  for  the  special  case  in  which  the  various  coef¬ 
ficients  satisfy  (29) ^  ^ *  They  are  the  same  as  those  noted  previously  in 
[12]  but  are  obtained  here  in  a  slightly  different  manner. 

Consider  now  a  typical  uniaxial  cyclic  loading  which,  for  the  ntl1 
cycle  in  the  s-e  plane,  passes  through  the  two  points  --  such  as  An  and  Dn 
in  Fig.  1  --  whose  (e,s)  coordinates  are  (nep,0)  and  (nep,0)  and  such  that 

!ne^|  >  |nep|  ,  min|ep|  =  |neP|  ,  max|eP|  =  |ne?|  ,  (33) 


where  the  notation  of  double  vertical  bars  signifies  absolute  value  of  the 


quantity  in  question,  and  where  min  and  max  refer  to  the  minimum  and  maximum 

values  of  plastic  strain  attained  during  the  nth  cycle.  Let  the  maximum  value 

of  compressive  and  tensile  stresses  attained  during  the  n**1  cycle  be  denoted 

respectively  by  "s^  and  and  similarly  during  the  n^  cycle  designate  the 

maximum  and  minimum  values  of  the  strain  respectively  by  "e^  and  "e^.  Then,  by 

(23)  the  maximum  and  minimum  values  of  the  plastic  strain  are  given  by 

1  n 


n  p  n  In 
*2  e2  "  E  S2 


n  p  n  In 

.  ei  *  ei  ”  F  si  •  (34) 

From  the  condition  f  «0  with  f  specified  by  (22),  the  uniaxial  stress  s 

on  the  yield  surface  is  found  to  be 
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f 


S  -  j  a(K)ep  t  V  J  <  , 


(35) 


where  the  plus  and  minus  correspond  to  tensile  and  compressive  values  of 
s  respectively. 

Next,  as  in  [4],  we  suppose  that  there  is  a  limiting  value  <s  of  the 
work-hardening  parameter  ic  as  t  ■*»,  i.e.. 


lim  <  =  k  ,  (36) 

t-*oo  S 

so  that  with  the  help  of  (18)  and  (25) 


lim  <  =  0 
t-*°° 

Then,  from  (35),  (13)  and  (34),  the  limiting  value  of  the  mean  stress  in  a 

* 

cycle  may  be  represented  as 


(37) 


_  — =  12  T  ("S1  '%>  *  r \  -V  -  F  '"=1 


lim  s 


(38) 


or  after  a  slight  rearrangement  as 

3Ea 

Smean  =  el  +Ile2)  =  ( 

t-*°°  s  n-*50 


3Ea 


-)e- 


4E+3a  mean 
s 


(39) 


where  for  cycling  between  fixed  strain  limits  (for  which  the  superscript 
n  may  be  suppressed)  we  have  defined 


e 


mean 


1  .  .  n  , 

j  lim(  e:  +  e2) 

n-»> 


(40) 


to  be  the  mean  value  of  the  strain  limits.  Moreover,  in  view  of  (18)  and 
(36),  from  (26),  (27),  (15)  and  (13)  we  obtain 


lim  $  = 

t-t» 


a 


s 

as«-4p 


<  1 


(41) 


Hence,  if  ag  >0,  then  by  (9)  and  (41)  the  material  saturates  to  the  value 


It  should  be  noted  that  for  repeated  cyclic  loading  the  number  of  cycles 
is  a  monotonically  increasing  function  of  time  so  that  in  the  present 
development  the  limits  at  t  ■*«  and  as  n  ■*<*>  can  be  interpreted  to  be 
equivalent . 
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given  by  the  right-hand  side  of  (41).  Correspondingly,  with  ag  =0,  the 
right-hand  side  of  (41)  vanishes  and  by  (9)^  the  material  saturates  to 
perfectly  plastic  behavior.  Since  eP/s  *dep/ds,  with  the  help  of  (23) 


and  (24)  we  obtain 


i  (s  -f  ae P) 


a(s  4aep)  ♦  [!  +(-S_l)(s  -j  aep)ep] (Bs  4  nep) 


K  -< 

0  s 


and  the  slope  of  the  stress-strain  curve  during  loading  is  given  by 

ds  .de.-l  A  dep,-l 
de  "  ldsJ  =  lE  +  dsJ 


At  initial  yield  this  has  the  value 


(— )  =  f—  + 

'•deJ  p  _  lE  a 

e^=0,tc=K  o 

o 

In  the  limit  as  t  it  follows  from  (43)  and  (42)  that* 


0  if  a  =  0  , 

s 


.  3Ea  9yklC 

de  *  4E+3ag  =  3k  lf  as  >  0  ’ 


where  in  writing  the  second  of  (45b)  use  has  been  made  of  (41)  and  (9a)  and 
the  elastic  constants  u  and  k  were  introduced  earlier  in  (11). 

To  demonstrate  the  predictive  capability  of  the  basic  theory  for  the 
type  of  response  described  under  1)  and  2)  of  Section  1,  we  consider  now  a 
uniaxial  elastic-plastic  state  specified  by 


,  eP  *  ep  t  0  ,  k  *  k 


and  suppose  that  the  material  is  unloaded  from  this  state  and  the  load  is 


reversed  until  it  again  yields  in  the  reverse  direction  at  s  =s. 


For  fixed 


values  of  ep  and  k,  s,^  and  s.jj  correspond  to  the  two  roots  of  (22)  given 


The  condition  (4Sa)  also  follows  from  the  first  of  (45b)  with  <*s  *0. 
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by  (35).  It  follows  from  the  sum  and  difference  of  these  roots  that 


*  1  *  *  2  ~  *  2  Sfll  *SP) 

K  =  F  Cs(l)  ' S (2) ^  *  a(<  }  =  I  ep*  •  ( 

A  value  of  <  follows  from  f  «0  at  initiation  of  yield  and  x  is  calculated 
o  1  s 

from  (47)  and  the  experimental  data  from  stress  reversals  at  saturation. 

It  follows  from  (39)  and  (45b) ^  that 

lim  s 

,  _  mean 

, .  ds  t-“> 

lim  j —  =  -  ,  (■ 

de  e 

t-*°°  mean 

so  that  a  value  for  the  quantity  a  may  be  determined  either  from  (45b)  and 

5  1 

the  measured  slope  of  the  stress-strain  curve  at  saturation,  or  from  (39) 

and  the  limiting  value  of  the  mean  stress  in  strain  cycling  between  fixed 

strain  limits.  An  elastic-plastic  state  between  the  initial  yield  and 

* 

saturation  is  specified  by  (46) .  After  evaluating  (15)  at  <  =<  and 
equating  the  result  to  (47),,,  there  results 


,  S,,,  +S.,.  X  -X  X  -X 

V  ^1-^'**^)  . 

eH  x  -x  x  -x 


where  <  is  calculated  from  (47) Next,  we  introduce  (49)  into  (44)  to 
determine  8  from  the  expression  for  the  slope  of  the  stress-strain  curve  at 
initial  yield.  Since  n  does  not  occur  in  the  yield  function  and  does  not 
appear  in  either  of  the  limiting  slopes  (44)  or  (45a, b),  its  value  must  be 
determined  by  using  the  general  expression  resulting  from  (42)  and  (43). 
Hence,  we  need  to  measure  the  slope  of  the  stress-strain  curve  at  several 
points  during  loading  between  the  initial  yield  and  saturation.  The  value 
of  x  at  such  a  loading  point  may  be  determined  from  (22)  and  (13),  or  from 
(47)  if  it  is  at  the  point  of  a  stress  reversal.  This  value  is  then  used  in 
(42)  and  (43)  to  estimate  an  appropriate  value  for  n. 
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The  above  procedure  is  suggested  as  one  way  of  determining  the  material 


constants,  and  is  by  no  means  unique. 

Fbr  purposes  of  illustration,  we  now  specify  the  material  constants 

for  a  material  similar  to  the  2024-T351  aluminum  alloy  considered  in  [4]. 

In  particular,  we  assume  the  same  values  for  E,xo  and  <s  as  those  given  in 

[4,  Eqs.  (92)  and  (93)].  In  addition,  we  specify  as  to  be  zero  so  that  the 

mean  stress  will  progressively  relax  and  tend  to  zero  during  strain  cycling. 

Now  instead  of  using  the  procedure  outlined  in  the  preceding  paragraph,  the 

coefficients  aQ,B  and  n  were  adjusted  until  the  predicted  stress-strain  curve 

exhibits  ratcheting  behavior  under  a  loading  program  of  stress  cycling  between 

* 

unsymmetric  stress  limits  .  In  this  way,  one  is  able  to  choose  the  material 
constants  as 
a 

-f«0.14  ,  as  =  0  ,  |=0.03  ,  4=-2.1*10"3  , 

<  <  . 

-y  *  1.8  xio'5  ,  —sS.Oxio'*  ,  (50) 

E  EZ 

E  =  69  GPa  . 

Stress-strain  response  curves  were  calculated  by  explicit  numerical  integra¬ 
tion  of  (24)  using  the  values  (SO)  for  the  material  coefficients.  The 
calculated  uniaxial  material  response  to  stress  cycling  between  fixed 
unsymmetric  stress  limits  is  shown  in  Fig.  1,  and  the  calculated  response  to 
strain  cycling  between  two  sets  of  fixed  strain  limits  is  shown  in  Fig.  2. 

From  (IS),  (26),  (35),  (27)  and  the  values  in  (50)  ,  the  value  of  4>  is 
calculated  to  be  0.10  at  initial  yield  and  it  tends  to  zero  at  saturation. 

Finally,  a  brief  comparison  is  made  of  the  prediction  of  the  theory  with 
experimental  results  for  304  stainless  steel  reported  by  Pugh  et  al.  {13]. 

This  manner  of  determining  the  coefficients  a0,8  and  n  is  analogous  to  the 
procedure  used  in  (4]  for  obtaining  different  coefficients. 

Also  assuming  a  value  v  »0.3  for  Poisson's  ratio,  with  ( 50)  the  shear 
modulus  p  is  calculated  to  be  26.5  GPa. 
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r 


Employing  the  procedure  described  earlier  fsee  the  paragraph  following 
(45)],  the  material  constants  were  determined  to  be 


ct  ,  a  .  „  - 

-j-  =1.4xio‘“  ,  =  1.7  xio'14  ,  |=3.2xio  , 

_  <  -  <  , 

*4  =  -1.2  xlO_,S  ,  -y=S.4xlO  ,  -4  =2.7x10  ,  (51) 

E  E  E 

E  =  123  GPa  . 

The  caluclated  reponse  is  shown  in  Fig.  3  along  with  the  experimental  curves 
for  the  first  three  cycles  from  [13],  as  well  as  the  results  calculated 
previously  by  Caulk  and  Naghdi  [4].  The  only  noteworthy  difference  between 
the  present  calculations  and  those  in  [4]  is  that  the  stress  level  at 
saturation  matches  the  experimental  data  even  more  closely  in  the  present 

* 

calculations.  From  (15),  (26),  (35)  and  (27),  (9),  (41)  and  the  values  in 
(51),  the  value  of  <t  is  calculated  to  be  0.029  at  initial  yield,  and  it  tends 
to  the  value  K,  =0.0053  at  saturation. 
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* 

As  in  [5]  assuming  a  value  v  «0.3  for  Poisson's  ratio,  u  is  calculated 
to  be  47.3  GPa. 
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Caption?  for  Figures 

Fig.  1:  Calculated  uniaxial  stress- strain  response  in  s-e  plane  to 

cyclic  loading  between  fixed  values  of  stress  in  tension  and 
compression  (not  necessarily  equal  in  magnitude)  resulting  in 
ratcheting  of  strain.  The  material  constants  used  in  the 
calculation  are  given  by  Eq.  (5(1). 

Fig.  2:  Calculated  uniaxial  stress-strain  response  in  s-e  plane  showing 
progressive  relaxation  to  zero  of  the  mean  stress  caused  by 
strain  cycling  between  any  two  fixed  values  of  strain.  The 
material  constants  used  in  the  calculation  are  given  by 
Eq.  (50). 

Fig.  3:  Comparison  of  theoretical  cyclic  stress-strain  behavior  for 

304  stainless  steel  with  the  experimental  data  of  reference  [13]. 

The  theoretical  stress-strain  curves  ( - ■)  are  calculated  using 

the  material  constants  in  (51) ;  the  theoretical  curves  (——•—) 
are  those  calculated  previously  in  [4];  and  comparison  with  the 

experimental  data  ( - )  is  shown  for  the  first  three  strain 

cycles  only,  since  curves  for  additional  cycles  would  crowd  the 
f i gure . 
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